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We study the low energy properties of the one-orbital quantum double-exchange model by using 
the slave fermion formulation. We construct a mean-field theory which gives a simple explanation 
for the magnetic and thermodynamic properties of the ferromagnetic metallic phase in manganites 
at low energy. The resulting electron spectral function and tunneling density of states show an 
incoherent asymmetric peak with weak temperature dependence, in addition to a quasiparticle 
peak. We also show that the gauge fluctuations in the ferromagnetic metallic phase are completely 
screened due to the Anderson-Higgs mechanism. Therefore, the mean-field state is robust against 
gauge fluctuations and exhibits spin-charge separation at low energy. 
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I. INTRODUCTION 

The phenomenon of colossal magnetoresistance (CMR) 
is known to occur in various doped perovskite manganese 
oxides with the chemical formula Rei_a:Aa:Mn03, where 
Re is the rare earth such as La or Nd and A is a divalent 
alkali such as Sr or Ca.— More recent studies revealed a 
complex phase diagram and very rich physicsi^ Therefore, 
it is of great theoretical interest to find a proper minimal 
model for these CMR systems, which can account for 
the important common features, such as the transport 
and magnetic properties, shared by all these materials, 
while leaving out many non-universal peculiarities due to 
crystal environment and atomic structure of individual 
compounds. 

One of the universally recognized common feature of 
these CMR compounds is the sizable ferromagnetic (FM) 
Hund's rule coupling Jh between the core spins and those 
of the eg electrons. This gives the double-exchange (DE) 
interactionSii^'^'^ii which is believed to be a fundamen- 
tal mechanism in explaining many of the interesting fea- 
tures of the CMR compounds. From this point of view, a 
usual starting point to analyze the properties of the CMR 
material is the standard DE Hamiltonian, supplemented 
with other terms, such as the Hubbard repulsion U, the 
antiferromagnetic (AF) superexchange interactions be- 
tween core spins Jaf, etc. For the case of CMR mangan- 
ites, band theory calculations^ suggest the typical values 
of the hopping amplitude t ^ 0.3 — 0.5 cV, Jh ^ 2.5 cV, 
and f7 - 6 - 8 eV. 

Motivated by experimental findings, there has been 
significant progress in the study of the DE Hamilto- 
nian, notably applying the Schwinger boson^iiS or 1/Sc 
expansion.— {Sc — 3/2 is the core spin.) Most 
of these approaches lead to a simple Fermi liquid (a 
doped band insulator) picture for the ground state of 
the FM metallic phase, which cannot fully explain the 
results of experiments, especially the optical conductiv- 



ity at low energyi^iii and the angle-resolved photoemis- 
sion (ARPES) measurements J^ii^iii Recently, Golosov 
tried to address parts of the discrepancy by including 
the on-site Coulomb repulsion through the Hartree-Fock 
approximationii^ Although some interesting results were 
found, the Hartree-Fock approximation is still based on 
the Fermi-liquid picture. To capture the low-energy 
physics for large values of U and Jh properly, a slave 
fermion approach has been proposed i^°'^^ A recent work 
by Hu2S along this direction has shown that the massless 
fluctuations of the longitudinal part of the gauge fields 
arising from the slave-fermion approach indeed dramat- 
ically change the behavior of the spectral function of Cg 
electrons at low energy. 

The purpose of the present article is to give a simple 
mean-field description on the low temperature proper- 
ties of the DE system in the large U, Jh limit. (For the 
implementation of the large U, Jh limit, see Sec. HIl) Or- 
bital fluctuations, Jahn- Teller effect, and nanoscale phase 
separation, though very interesting, make analysis diffi- 
cult. Therefore, we shall focus our attention on the region 
where the DE mechanism is the dominant factor, namely, 
the region at the hole concentration 0.2 < a; < 0.5 and 
away from the critical temperature Tc- We propose a 
mean-field theory based on the slave fermion scheme, 
which exhibits spin-charge separation at low energy, and 
use it to calculate various low-energy properties in the 
FM metallic phase. Among them, the most important re- 
sult is that the quasiparticle peak in the electron spectral 
function is reduced at low temperatures and the spectral 
weight is transferred to an asymmetric broad peak away 
from the Fermi surface, which is a natural consequence of 
the spin-charge separation in our theory. The other phys- 
ical properties in the FM metallic phase we have studied 
are as follows: (i) The magnon dispersion at the mean- 
field level is similar to that of a simple cubic Heisenberg 
ferromagnet, which has been verified by inelastic neu- 
tron scattering measurements, (ii) The magnitudes of 
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the spin stiffness at T = and the coefficients in the low- 
temperature specific heat, obtained from the mean- field 
theory, are consistent with experimental data, (iii) The 
structure of the optical conductivity at low energy is of 
the form similar to that observed by experiments, and 
the Drude weight is reduced. 

We now briefly outline the structure of our paper. In 
Sec. ini we will introduce the quantum double exchange 
(QDE) model and develop a slave fermion mean-field the- 
ory. The results of the mean-field theory is presented in 
Sec. mil We will show the stability of our mean-field 
state against the gauge fluctuations in Sec. [IVl The last 
section is devoted to a conclusive discussion. 



II. THE MODEL AND THE SLAVE FERMION 
MEAN-FIELD THEORY 

We shall start with the quantum double-exchange 
(QDE) model described by the following Hamiltonian>i^ 



H.c.) - JH^5?-Si , (1) 



where the flrst and second terms describe electronic hop- 
ping and DE couplings, respectively. Here, S'^ denotes 
the core spin, s is the spin operator of Cg electrons, m is a 
unit vector connecting the nearest-neighbor sites around 
site i (for a simple cubic lattice, u — x, y, and z), and 
Cia is the annihilation operator of Cg electrons at site i 
with spin a. In Eq. ([T|), we have neglected the superex- 
change interactions between core spins. This is because 
Jaf ~ 5 — 10 meV in manganites, which is much smaller 
than t. Moreover, in the undoped compounds, a Jahn- 
Teller distortion lifts the orbital degeneracy of the Cg elec- 
trons with the energy scale Ejt ~ 1 — 1.6 eV. Thus, to 
study the low energy physics in the hole-doped region 
(x < 0.5), one may in a flrst approximation ignore the 
orbital degeneracy and apply the one-orbital model. To 
take into account the fact that the on-site Coulomb inter- 
actions between Cg electrons is the largest energy scale in 
manganites, we impose the no-double-occupancy (NDO) 
condition on the Cg electron operators: 



0,1 



(2) 



Therefore, Ci„ and ct^ cannot obey the canonical anti- 
commutation relations, and Eq. ([T]) is valid for energies 
much lower than U . 

Because of the NDO condition [Eq. the system 

becomes a strongly correlated one. One popular way to 
solve this kind of problems is to introduce a pair of slave 
flelds to rewrite the c operator. Here we adopt the slave- 
fermion representation: 



where fi and fj satisfy the canonical anticommutation 

relations, and bia and satisfy the canonical commuta- 
tion relations 1^ In this way, the charge and spin degrees 
of freedom of Cg electrons are represented by / (holon) 
and b (spinon) fields, respectively. In terms of fi and bi^, 
the NDO condition becomes an identity 



Y.blbi. + flfi ^ 1 



(4) 



In manganites, one may further simplify Eq. ([T]). Since 
Jh ^ lOi, at the energy scale much lower than Jh, it suf- 
fices to consider the Hilbert space in which the Cg elec- 
tron spin is parallel to the core spin. One may intro- 
duce another Schwinger bosons dia and to describe 
the total spin. Within this subspace, it can be shown 
that the b field is associated with the d field through the 
relationiSiii^ 



(5) 



/2S 



with 5 = 2. Collecting the above results, in the limit 
jR/t —> +00, Eq. ^ is reduced to the one 



Hde = XI + l^-c-) 



(6) 



and the NDO condition becomes 



(7) 



We shall take Eq. (I6|) as a minimal model to describe the 
FM metallic phase in manganites.— 

To proceed, we turn into the path-integral formalism. 
The partition function of the QDE model in the large Jh 
limit can be written as 



Z = jD[p]D[f]D[di]D[d,]D[X]expl^-jyTj2L 
where 

^ = I]rfL(9r + Ai)di, + /|(9,+zAi-Mo)/i 



(8) 



(3) 



fiQ is the chemical potential of holes, and Aj is the La- 
grangian multiplier to impose the NDO condition. To 
facilitate the mean-field analysis, one may perform the 
Hubbard-Stratonovich transformation to decouple the 
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hopping term, and the Lagrangian becomes 

L = ^ dlai^T + iXi)dia + fl{dr + iXi - l-Lo)fi 



2S X-/ ("*''^+"'*-^*^+"'^'' ^ 

u,cr 

I E + H.c.) - 2i5Ai, (9) 



We notice that there is a U(l) gauge structure in the 
slave-fermion scheme, which is reflected in the invariance 
of L (up to a total derivative in t) under the U(l) gauge 
transformation: 



f. ^ f.p-™i rl. ^ fj. p- 



Ai ^ Ai + 9rw;i 



(10) 



Motivated by the DE mechanism, which allows the co- 
herent hopping of charge carriers in a FM background, 
we consider the following mean-field ansatz: 



Xi+u.i = X , Vi+uA = V , 



(11) 



where x, rj, and A are real. Inserting Eq. into Eq. 
([9|) , the mean-field Lagrangian can be written as 



H.c, 



25 

U 



where z = 6 is the coordination number. Now Lmf be- 
comes quadratic in / and d^^, and one may integrate 
them out to obtain the mean-field grand potential. The 
mean-field equations are obtained through minimizing 
the mean-field grand potential with respect to the mean- 
field parameters, yielding 

^Y.^Bie,k + A) = (12) 

k 



zt 



N ^ 

k 



ebknB{f-bk + A) = —XV 



S 
zt 
25 



XV 



(13) 
(14) 



where np{x) = (e^^ + 1) ^ is the Fcrmi-Dirac distri- 



bution, ng^x) 



1) ^ is the Bose-Einstein dis- 



tribution, is the number of lattice points, eff. 



Y.U cos {k-u), and e^fc 



f E.cos(fe.M). The 



chemical potential of holons, /x^ = /.tp — A, is associated 
with the hole concentration x through by the equation 



]^Y.''F{efk 



M/) 



(15) 



Equations p2|) — p5|) are the mean-field equations we 
want to solve. 

Before solving these mean-field equations numerically, 
we may look into the physics revealed by them. Equation 
(|12p is nothing but the number equation of free bosons 
with the average number of particles per site S — x/2 
and the chemical potential —A. In three dimensions, 
there exists a critical temperature Tc > so that A = 
Ao for T < and A > Aq for T > where Aq = 
T]zt/{2S) is the bottom of the spinon band ebk- When 
A = Aq, the Bose-Einstein condensation (BEC) of the d 
bosons occurs, which corresponds to the FM phase. On 
the other hand, it is the paramagnetic (PM) phase for 
A > Aq. In this scenario, the PM to FM phase transition 
in manganites is associated with the BEC of the d bosons, 
and Tc corresponds to the Curie temperature. 

With the above understanding in mind, we can see that 
deep inside the FM phase, without loss of generality, one 
may choose the direction of magnetization to be the z 
axis and set {di-\) = \/2S — x and (di^) = in terms of 
a proper spin SU(2) rotation. (This parametrization for 
d bosons is consistent with the experiment which shows 
the complete spin polarization of conduction electrons in 
manganiteS;--) Further, for T ^ Tc, one may neglect 
the amplitude fluctuations of As for the phase fluc- 
tuations of di-\ , it can be absorbed into dii and fi by 
choosing a particular gauge or performing a proper U(l) 
gauge transformation with the help of the U(l) gauge in- 
variance of L. Thus, to obtain the low-energy effective 
Hamiltonian in the FM phase, it suffices to set 



1^^ 



V2S'- 



(16) 



in Lraf- On account of the condensation of the d bosons, 
it turns out that the gauge fluctuations acquire a finite 
energy gap through the Anderson-Higgs mechanism, so 
that the mean-field state is stable against the gauge fluc- 
tuations. As a result, the physics in the FM phase at 
energies much lower than the gap of gauge bosons, Eg^ 
can be described by the following effective Hamiltonian 
in the grand canonical ensemble: 



Hfm = E \_^^^ ^ t^f)fifk + '^kbihk 



(17) 



where h = ^Hk and hi = ^Y.k ^^'^^bk- (For 

the details of the derivation, see Sec. IIVI ) The / field 
describes the spinless charged excitations — holons, with 
the dispersion relation 



(18) 
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and we shall see that the b field describes the FM spin 
waves — magnons, with the dispersion relation 



rjt 
2 



3 — ^ cos {k ■ u) 



(19) 



In Eq. ((TT)) . we have neglected the interactions between 
holons and magnons, which are irrelevant operators in 
the sense of the renormalization group (RG). Equation 
(|17|) indicates that there is spin-charge separation at low 
energy in the FM metallic phase. (We will come back to 
this point in Sec. II VI ) Further, we shall sec later that 
the bandwidth of holons is of the order of t, while that 
of magnons is of the order of O.lt. 



III. RESULTS IN THE FERROMAGNETIC 
METALLIC PHASE 

A. Magnetic properties 

a. Magnon dispersion We shall first employ our 
mean-field theory to calculate the transverse spin-spin 
correlation function, which is defined as 



(20) 



where Sf — Sf ± iS^ and x = i — j . In the large Jh 
limit. S± can be expressed by the d bosons 



ST 



(21) 



The Fourier transform of S±{t, x), denoted by S±{uj, q), 
can be obtained from the corresponding Matsubara func- 
tion through analytical continuation. In the FM phase, 
S±{lo, q) can be related to the two-point correlation func- 
tion of the h field within the mean-field theory: 



S^{LO,q)^Xd{T)S'^{Lo,q) 



(22) 



where Xd{T) = {d\-^di-^) is the average number of d bosons 
in the condensate per site at temperature T and S'^{uj, q) 
is the Fourier transform of the retarded Green function 
of the b field. In view of Eq. , one may identify the 
excitations corresponding to the b field as the magnons. 

We notice that the form of the magnon dispersion we 
obtained [Eq. ((T9)) ] is identical to that given by the FM 
Heisenberg model on a simple cubic lattice with nearest- 
neighbor interactions only. This result is similar to that 
predicted by previous studies on the DE modeli^i^'^ii 
However, within the present framework, the Heisenberg- 
type behavior for the magnon dispersion is a mean-field 
result. By taking into account the so far ignored ir- 
relevant (in the sense of the RG) interactions between 
magnons, a deviation from the Heisenberg spectrum is 
expected. On the general ground of the RG, we expect 
that the deviation from the mean-field result (the Heisen- 
berg spectrum) will become more noticeable away from 



the zone center. Furthermore, the small values of the 
spin stiffness at T = in the regions a; — *■ and x — > 1 
(see Fig. [1]) suggest that the mean-field results may re- 
ceive considerable corrections in these regions. On the 
other hand, for the physically interested doping regime 
0.2 < a; < 0.5, the mean-field results should be robust. 
Recent works on the DE model, such as the spin-wave 
theory based on the large Sc expansio n^^i^^ and an exact 
calculation on a finite ring2i in the limit Jh /t ^ 1, have 
revealed a deviation of the magnon dispersion from the 
Heisenberg spectrum in the DE model. The deviation is 
prominent at large momenta^i^ and the overall devia- 
tion is very small in the doping range 0.2 < x < 0.6^ 
Our theory is consistent with these results. 

A Heisenberg-type behavior for the magnon dispersion 
may provide a reasonably accurate picture for manganese 
oxides with large values of However, recent exper- 
iments indicate deviations from this canonical behavior 
in compounds with lower values of Tc- In particular, the 
softening of the mangon dispersion near the boundary 
of the Brillouin zone is observed Such an issue clearly 
depends on the details of the short-distance physics, and 
is beyond the scope of the present work. A systematic 
calculation which incorporates the magnon-magnon and 
magnon-holon interactions, such as those given by Eq. 
p7)) . may be helpfuli^ However, for a complete compar- 
ison with experimental data, additional ingredients, such 
as orbital fluctuations and orbital-lattice couplings, may 
also need to be taken into account 

A finite damping rate of magnons at T = in the 
DE model has been pointed out in Refs. [HI [13, and[28l. 
To address this problem within the slave-fermion theory, 
we must go beyond the mean-field results. In fact, the 
magnon acquires a finite lifetime, even at T — 0, at the 
two-loop order of the magnon self-energy diagrams by 
taking into account the interactions between holons and 
magnons. [The leading terms of these interactions are 
given by Eq. (j47p .] The physical origin of a nonvanish- 
ing T = damping rate of magnons can be easily un- 
derstood based on our spin-charge separated mean-field 
state: Because the magnon dispersion, which is propor- 
tional to near the zone center, is definitely immersed 
into the "particle-hole" continuum of holons in three di- 
mensions, the magnon can decay by exciting a single 
"particle-hole" pair and another magnon, which is similar 
to the physics of Landau damping in the Fermi liquidsi^ 
Since the interactions between holons and magnons are 
irrelevant operators in the sense of the RG, we expect 
that the damping rate of the magnon near the zone cen- 
ter (fc = 0) is, at least, proportional to k" with a > 4 
in three dimensions by the power counting usually cm- 
ployed in the momentum-space RG. (Additional on-shell 
constraints may further increase the value of a.) The 
spin- wave theory in the Jh /t +oo limit predicts that 
the damping rate of the magnon near the zone center is 
proportional to k^ in three dimensions fiiiiS which implies 
that it is indeed due to the irrelevant interactions. A full 
calculation by incorporating these irrelevant operators is 
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FIG. 1: (Color online) The doping dependence of the spin 
stifltness at T = 0. Here we set the lattice constant 1 = 1. 



FIG. 2: (Color online) (a) Tc/t v.s. x; (b) fif/t v.s. x at 
T = Tc; (c) »7 v.s. a; at T = Tc, (d) x v.s. a; at T = Tc 



beyond the scope of this paper, and their effects within 
the present framework are reserved for further studies. 

b. Doping dependence of the spin stiffness at T — 
From Eq. ([TO]) , the magnon dispersion in the long 
wavelength limit fc ^ is given by ujk ~ DqU^ where 



ritP 
'2S 



(23) 



the present mean-field theory, and it is 



M{T) = Mo 



T 



C(3/2) 
Afo V 47rL»o 



3/2 



(24) 



which is identical to that predicted by the simple cubic 
Heisenberg ferromagnet, where Mq = (2 — x/2)l^^ is the 
magnetization at T = 0. 



is the spin stiffness at T = and I is the lattice constant. 
Thus, the doping dependence of can be extracted from 
that of ?/. The result is shown in Fig. [1] 

Our theory predicts that Dq is symmetric in x with 
respect to the quarter-filling x — 0.5. This result is sim- 
ilar to that predicted by the spin-wave theory of the DE 
model in the limits Jn/t +oo and Sc +oor^ ex- 
cept that the magnitude of Do{x) around x = 0.5 we 
obtained is smaller. The values of Dq we obtained are 
Do/{tP) = 0.0144,0.018,0.0202 for x = 0.2,0.3,0.4. Us- 
ing the relation Dq = JSeffl^, where Seff = 2 — x/2 
is the average spin for each site and J denotes the ef- 
fective exchange coupling between spins, one may get 
J = 0.0076i,0.0097i,0.0112t for x = 0.2,0.3,0.4, which 
correspond to J = 2.28, 2.91, 3.36 meV if we use t = 0.3 
eV. J w 1.9,2.4 meV found in Lao.8Sro.2Mn03^ and 
Lao.rPbo.aMnOar^ respectively. We see that the values 
of J we obtained are consistent with those extracted from 
experimental data. To sum up, for the magnon disper- 
sion relation, the difference between our mean-field the- 
ory and the spin-wave theory of the DE model in the 
doping range 0.2 < x < 0.5 mainly lies at the doping 
dependence of the spin stiffness at T = and its magni- 
tude. 

Finally, we mention that the low temperature magne- 
tization in the FM phase can be easily calculated within 



B. Estimation of Tc 

Next, we would like to estimate Tc within our mean- 
field theory, which can be obtained from the mean-field 
equations [Eqs. (fT2|) — ([T5|) ] numerically by setting 
A = Aq. The mmicrical results are shown in Fig. [2l A 
few points about our results should be discussed. First 
of all, Tc is asymmetric in x with respect to the point 
X = 0.5 though T][Tc) is still symmetric in x with respect 
to a; = 0.5. This is in contrast to Dq. This asymmetry 
suggests that the FM phase is more robust at a; < 0.5 by 
including strong on-site Coulomb repulsions and Hund's 
rule couplings properly. Next, the previous work based 
on the dynamical mean-field theory gives an estimate of 
Tc for Jn/t > 1: Tc/t = 0.146 at a; = 0.3,1 Our mean- 
field theory predicts Tc/t ~ 0.1454 at x — 0.3, which 
is quite close to the value obtained from the dynami- 
cal mean-field theory. Finally, the maximum of Tc/t is 
reached at a; = 0.458 with the value Tc/t = 0.1575, which 
corresponds to Tc ~ 548 K for t ~ 0.3 eV. This value 
is about 2 times of that obtained by experiments. As 
usual, the fluctuations will reduce the mean-fleld value. 
Nevertheless, the difference between the mean-field result 
and experimental data may still not be explained even if 
we include the fluctuations. This may not be surprising 
since, after all, Tc is a non-universal quantity and the 
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FIG. 3: (Color online) 7 = c^-uc/T, in unit of mj/mol K'^, in 
the doping range 0.2 < x < 0.5. 



FIG. 4: (Color online) B3/2 = Cmag/T'^^'^, in unit of mJ/mol 
K^''^, in the doping range 0.2 < a: < 0.5. 



above result simply indicates that to have a good esti- 
mation of Tc a few ingredients that are ignored in the 
QDE model, such as the AF superexchange interactions 
between core spins, orbital fluctuations, and electron- 
phonon interactions, have to be included. In fact, ex- 
periments on the oxygen-isotope substitution show that 
phonons arc important in the determination of Tc.— 



C. The low-temperature specific heat 

Low-temperature (T < 10 K) heat-capacity measure- 
ments provide information regarding the bulk properties 
of solids. For a magnetic solid, the low-temperature spe- 
cific heat is composed of numerous contributions and it 
is typically given by 



(25) 



Here Ceiec is the electronic contribution, which takes the 
form Ceiec = iT , Chyp arises from the hyperfine field of the 
nuclear moment, ciat is the contribution of phonons, and 
is the contribution from the magnetic spin waves. 



^7nag 



which is usually estimated as the form J2n ^nT". 

In our theory, Ceiec is primarily given by the holon 
sector due to the spin-charge separation at low energy, 
yielding 



7 



away from the van Hove singularities, where fc^ is the 
Boltzmann constant, and giep) is the density of states 
(DOS) of holons at the Fermi energy. 

Figure [3] shows the values of 7 in the doping range 
0.2 < a; < 0.5 with t = 0.3 eV (solid line) and t = 0.2 
eV (dashed line). The magnitudes of 7 we obtained are 



consistent with the experimental data. Furthermore, the 
mean-field theory predicts a very weak doping depen- 
dence of 7, which is approximately given by 



7 oc 



1-x/A 



(26) 



Such a weak doping dependence of 7 may not contra- 
dict to the experimental data at x = 0.2 and 0.3, where 
the measured values of 7 (« 3.3 mJ/mol K^) appear to 
be nearly a constant in the FM metallic phasci^ More 
detailed experiments is warranted to resolve this issue. 



As for 



the magnon contribution to the low- 



temperature specific heat is similar to that in the simple 
cubic Heisenberg ferromagnet and takes the form 



ith 



-'mag 



-B3/2 — 



15C(5/2)4/" 



The values of i?3/2 in the doping range 0.2 < x < 0.5 are 
shown in Fig. |4]with t = 0.3 eV (solid line) and t ~ 0.4 
eV (dotted line). The values of B3/2 we obtained are close 
to the experimental data^ (-B3/2 = 1.1802 mJ/K^/^ mol 
for X = 0.2 and B3/2 = 0.95672 mJ/K^/^ mol for x = 0.3) 
if we use t = 0.4 eV. 

To sum up, the consistency of the magnitudes of 7 
and -83/2 in the doping range 0.2 < x < 0.5 evaluated 
in terms of the mean-field theory suggests that the low- 
temperature thermodynamics of the FM metallic phase 
can be well-described by the one-orbital QDE model, and 
this result can be viewed as an indirect evidence of spin- 
charge separation at low energy in the FM metallic phase. 
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spectra are composed mostly of the incoherent part and 
hghtly of the Drude response (about 20% - 30% in frac- 
tion), (ii) The Drude part is discernible below 0.04 eV, 
but with an anomalously small spectral weight. For ex- 
ample, the value of DW is as small as 0.012 even for 
the lowest temperature spectra for both x = 0.175 and 
0.3. Further optical studies on Lao.rCao.sMnOai^ con- 
firm that, in the FM metallic phase at very low tem- 
peratures, the low-energy optical conductivity spectra 
{lo < 0.5 eV) show two types of absorption features: a 
sharp Drude peak with little weight (about 33% in frac- 
tion or DW ^ 0.02) superimposed to the broad incoher- 
ent absorption band. Our analysis indicates that in the 
FM metallic phase the optical conductivity spectra of the 
QDE model at low energy {u; ^ Jh) indeed consists of 
two parts: 



FIG. 5: (Color online) DW as a function of x. 



(71 (w) CTc(cj) -I- (Jinci^^) , 



D. Optical conductivity at low energy 

One of the consequences of the spin-charge separation 
in our theory is that the optical conductivity at low fre- 
quency (w ^ Eg) is mainly contributed by free holons, 
where Eg denotes the gap of gauge fluctuations. Espe- 
cially, the Drude weight is completely determined by free 
holons, which is given by 

= (27) 

where K is the average kinetic energy of holons per site 
and I w 3.9 A is the lattice constant. Following exper- 
iments, one may measure the Drude weight in unit of 
7re^/(2meZ^), denoted by DW, where rrie is the electron 
mass. Within our mean-field theory, DW at T = is 

W.5!I^(l_£), (28) 

where h is the Planck constant. The doping dependence 
for DW, given by Eq. (|28p . is shown in Fig. O with 
t = 0.3 eV (sohd hne) and t = 0.2 eV (dotted line). 
From it; we see that, in the doping range 0.2 < a; < 0.5, 
DW ~ 0.19 - 0.26 for t = 0.3 eV, and DW - 0.13 - 
0.18 for t = 0.2 eV. Since the Drude peak at T = 
does not carry 100% of the weight, this result implies 
that some spectral weight must be transferred to higher 
energies {lo > Eg) due to the optical conductivity sum 
rule, so that the Drude weight is suppressed compared 
with normal metals. 

Optical conductivity spectra have been investigated for 
single crystals of Lai-^jSr^jMnOa with < a; < 0.3fi^ 
The peculiar behaviors observed in the low-energy opti- 
cal spectra (w < 0.1 eV) in the FM metallic phase at low 
temperature, which cannot be explained by the simple 
Fcrmi-liquid picture, arc as follows: (i) The low-energy 



where adto) is the coherent (Drude) part, which is mainly 
contributed by holons, and aindio) represents the inco- 
herent part with most of its weight in the region to > Eg, 
which arises from the strong scattering between gauge 
fields and holons. Thus, according to our theory, the en- 
ergy scale below which the Drude part is discernible can 
be regarded as the lower bound oi Eg. A rough estimate 
of the value of Eg (Eg - 0.038 - 0.071 eV in the dop- 
ing range 0.15 < x < 0.5 for t ~ 0.3 eV) is consistent 
with the experiment. (See Sec. IIVH This picture is con- 
sistent with the observed optical spectra at low energy. 
However, the values of DW we obtained are larger than 
the experimental data by an order of magnitude. Since 
our mean-field theory is supposed to be accurate at low 
energy, the quantitative discrepancy between our results 
and the experimental data indicates that to explain the 
optical spectra of manganites, other ingredients beyond 
the one-orbital QDE model must be taken into account. 
Some possible scenarios, such as orbital fluctuations^ or 
electron-phonon conplingSf^i^i have been proposed. Ap- 
plying our method to these "extended" QDE models is 
an interesting problem. 



E. The electron spectral function and tunneling 
density of states 

The most important result in our mean-field theory is 
the spin-charge separation at low energy. A direct exami- 
nation of this phenomenon is to study the behavior of the 
electron spectral function at low temperatures, which can 
be measured by the ARPES experiments. The electron 
spectral function is defined as 

A{LU,q) = -2lm{g{iLUr,^uj,k)} . (29) 

Here Im{- • • } means the imaginary part of • • • , and g{K) 
with K — {iLOn, k) is the Fourier transform of the electron 
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Green function, which is defined by 

= -^{nfKr)dMdimjm) , m 

where X = {T,i — j). 

At the mean- field level, holons and magnons are decou- 
pled from one another and dii is treated as a c- number, 
in stead of an operator. As a result, the electron Green 
function can be approximated as 



(a) kp=[0.747i,0.74jt,0] (b) kp=[0.547i,0.547i,0.547i] 



GiX) 



+ ^{nmflir)}){ndiiir)dliO)}) 



2S 



G{-X) 



1 

25 



T"iT/ 
G{~X)S{X) 



where Xd{T) is determined by the equation 



(31) 



(32) 



and 



G{X) = -{T{h{T)f]m) , 

s{x) = -{T{b,{T)b]m) , 

arc the propagators of holons and magnons. respectively. 
In terms of Eq. ((3T|) . the electron spectral function is 
given hy^ 

A{oj, k) = Z{x, T) 5{uj + €k- f^f)+ a{uj, k) . (33) 
where Z{x,T) — 'KXd{T)/S and 

<? 

X 5{UJ + eq - flf - UJk+q) ■ 

Equation ([55]) indicates that the electron spectral func- 
tion at low energy consists of two parts: a sharp quasipar- 
ticle peak, following the holon dispersion, superimposed 
to a broad incoherent part given by a(a;, fe). The spec- 
tral weight of the quasiparticle peak is given by Z{x,T), 
which is associated with the condensate of the d bosons. 
At T < Tc, it is given by 



Z{x,T) w 2tt 



X 

4 



C(3/2) / TP 
4 K^ttDq 



3/2' 



(34) 



By taking into account impurity scattering and the in- 
teractions between holons and magnons, the (5- function 
peak will be broadened, and the width 1/r is given by 
1/r = maa;{l/ro, /cbT} where tq is the lifetime due to 
impurity scattering. The low-temperature behavior of 
a(uj,kp) at the hole concentration x — 0.4 is shown in 
Fig. M where kp « 0.747r(l, 1, 0) and 0.547r(l, 1, 1) in 
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FIG. 6: (Color online) The electron spectral function, with 
the subtraction of the quasiparticle peak, at the hole con- 
centration X = 0.4. kp denotes the Fermi momentum. The 
temperatures we consider are T/t = 1/174 (dashed line) and 
25/1088 (solid line), which correspond to T = 20,80 K, re- 
spectively, for t — 0.3 eV. 



(a) x=0.4 



(b) x=0.3 




FIG. 7: (Color online) The electron spectral function, with 
the subtraction of the quasiparticle peak, at the temperature 
T/t = 1/174 and hole concentration (a) x — 0.4 and (b) 
X = 0.3. Here k = ko x [n, tt, 0] with ko = 0.5, 0.6, 0.7, 0.8, 0.9 
from bottom to top. 



Fig. m (a) and (b), respectively. Figure [7] shows a(w,fc) 
along the direction from [0,0,0] to [tt, 7r,0], at the tem- 
perature T/t = 1/174 and the hole concentration (a) 
X = 0.4 and (b) x = 0.3. A common feature of a(w, k) 
is that it exhibits a broad asymmetric peak at low tem- 
perature with the width ~ t, and the peak position is 
at w ~ 0.5t and 0.21^ for x — 0.4 and 0.3, respectively. 
Furthermore, at low temperature, the peak position and 
its width are insensitive to the variations of tempera- 
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tures. This asymmetric peak implies the violation of the 
particle-hole symmetry. As raising the temperature, we 
can see a transfer of the spectral weight from the quasi- 
particle peak to the incoherent part. Such a trend per- 
sists until at the critical temperature where the coherent 
part disappears completely, i.e. Z{x,Tc) = 0. The pres- 
ence of such a broad asymmetric peak at low temperature 
reflects the composite nature of electrons at low energy 
in the large [/, Jh limit, i.e. a manifestation of the spin- 
charge separation at low energy, which can be observed 
in the photoemission experiment. 

A similar form for the electron spectral function [Eq. 
(j33p ] was also obtained by a mean-field approximation 
to the DE model in the limit Jn/t +00.— There, the 
FM phase also corresponds to the condensed phase of 
the Schwinger bosons. However, the doping dependence 
of the mean-field parameters in that approach is differ- 
ent from ours, especially the bandwidth of holons. This 
results in distinct doping dependence for various physi- 
cal quantities, such as the spin stiffness at T = and 
the Drude weight. Moreover, quadratic approximations 
for the dispersion relations of holons and magnons were 
used to calculate the electron spectral function, which 
leads to a weak logarithmic singularity in the incoher- 
ent part a(w, k) with the width being of the order of the 
maximum magnon energy. We do not see such a singu- 
larity, as can be seen in Figs. [n]and[71 Further, the width 
of the broad peak we obtained is much larger, which is 
consistent with the ARPES measurements. 

Another quantity which is related to A{uj, k) is the 
tunneling DOS: 

^H = ^E^(^'^) (35) 

k 

which can be measured by the scanning tunneling spec- 
troscopy (STS). Within the mean-field theory, N{uj) can 
be written as 

N{co) = ^^^^ nJ ^-)+ - Y a{uj,k) , 

(36) 

where N^ie) = ;^ I]fc '^['^ " Z^u ^os (fc • u)] and uj ^ 
{lu — fJ.f)/t. Figure [5] shows N{uj) as a function of Cu 
at two temperatures T/t = 1/174 and 25/1088 with the 
hole concentration (a) x = 0.4 and (b) x = 0.3. A few 
salient features about the low-temperature behavior of 
the tunneling DOS are as follows: (i) N{u}) consists of two 
parts: One is proportional to the DOS for a tight-binding 
model on a simple cubic lattice, which is contributed by 
the quasiparticlc peak of the electron spectral function, 
and the other exhibits an asymmetric peak, which re- 
sults from the incoherent part of the electron spectral 
function, (ii) The position of the asymmetric peak is at 
w — /i/ ~ xiT = 0)t, the renormalized hopping amplitude 
of holons. (iii) The values of N{uj) is insensitive to the 
variation of temperatures. Again, this structure of N{uj) 
is a natural consequence of spin-charge separation. 



(a) x=0.4 (b) x=0.3 




-5 5 -5 5 



CO CO 

FIG. 8: (Color online) The tunneling DOS at the tempera- 
tures T/t = 1/174 (solid line) and 25/1088 (-I-) with the hole 
concentration (a) x = 0.4 and (b) x — 0.3. We have set t = 1. 



Variable temperature STS studies on single crystals 
of Lao.rPbcaMnOs have been done in the temperature 
range 100 — 375 Ki^ Experimental data show that N^ep) 
grows rapidly below Tc and reaches a temperature inde- 
pendent value, where is the Fermi energy. Further- 
more, N(lu) near is flat in the scale of T below T^,. Fur- 
ther STS studies on single crystals of Lao.yCao^aMnOa^ 
confirm these results. Our prediction about the behavior 
of N{uj) near the Fermi energy is consistent with exper- 
imental results. Further STS studies at low tempera- 
tures must be performed to verify our mean-field theory. 
Especially, an observation of the asymmetric peak away 
from the Fermi energy in the tunneling DOS at low tem- 
peratures will be solid confirmation of the slave-fermion 
theory. 

IV. GAUGE FLUCTUATIONS AND THE 
STABILITY OF THE MEAN-FIELD STATE 

In order to study the role of gauge fluctuations and to 
derive the low energy effective theory, we can start with 
Eq. ([H]), and parameterize the Hubbard-Stratonovich 
fields as 

Xi+^,i - xe"'*^'" , m+^.i = 'ye-*- . (37) 
One may further decompose and (j)i_u by 

4^i,u -^i.it ^ ^i.u 7 4^i,u -^i.-u ^i,u ■ (38) 

Such a decomposition becomes transparent once we real- 
ize that under the U(l) gauge transformation, Ai,u trans- 
forms like a gauge field, while Bi^u is gauge invariant. 
Simple manipulations show that the Bi^u field acquires 
a finite energy gap. Furthermore, the amplitude fiuc- 
tuations of Xi+u.i and rji^^,! are also gapped. For the 
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energy below these mass gaps, one may integrate out the 
massive fields and take the continuum limit 

Ai,u ^ u ■ a{X) , Ai ^ X{X) , (39) 

where X = {t,x). After doing so, we arrive at the fol- 
lowing continuum Lagrangian within the effective-mass 
approximation 



C = 



J2 €{dr + iX)da + h\dr + i\ - ^ia)h - 2iMX 

a 



where M ~ Sl^^ . In view of Eq. (|40p . the Lagrangian 
multiplier A plays the role of the time component of the 
gauge fields. 

In the FM phase, without loss of generality, one may 
parameterize d„ in the following way: 



(41) 



on account of the BEC of the Schwinger bosons, where 
Ud = {2S — x)l~^ is the average density of d bosons, and 
p and (f> describe the amplitude and phase fluctuations 
of respectively. That is, we choose the direction of 
the order parameter (magnetization) to be the z-axis. 
Inserting Eq. (|4ip into Eq. (jlQ]) gives rise to 

^FM = h\dr + iX)b + h\dr + iX — ^.i^)h — in\ 



1 



2m/ 
1 



{V -ia)h\' 



2mb 



\{V~^a)b\^ 



nd 



ZiJib 

The U(l) gauge structure of Cfm now reads 

h ^ he-'"" , b 6e^™ , 

p , cj)^ (p-w , 
a ^ a — Vw , A — !• A -f drW . 



(42) 



With an eye on the gauge invariance of Cfm, one may 
choose the gauge w — <j). This amounts to redefining the 



h, b, a and A as follows: 



h = he-"'' , b = be-"!' , 
a = a — V0 , A = A -I- ^^(j 



(43) 



We notice that h, b, a, and A are all gauge invariant. In 
terms of Eq. (|43l) . Lfm becomes 



Cfm ~ b^dr + iX)b + h^dr + iX — po)h + i{p — n)X 



1 



■\{V -id,)hf' 



1 



2mb 



\iV-ia)b\' 



2m/ 

^ |Vpp + -^(nd-f p)a2. 



In the above, the higher order terms in p, i.e. the self- 
interactions of p, have been neglected. This is because 
they only generate irrelevant operators. Integrating out 
A gives rise to the constraint 

Pb + h^h + p = n , (44) 

which is nothing but the continuum version of the NDO 
condition [Eq. ([7])]. Using Eq. (jH]), one may further 
integrate out the p field, yielding 

Cfm = h^{dr- Po)h + ^\Vh\^ + b^drb + Do\Vb\^ 



2m 



/ 



1 -2 

2 '^aCl 



where • • • represents the interactions between h, b, and 
a, Dq = {2mb)~^ is the FM spin stiffness at T = 0, and 
nia = [2M — (1 — mb/mf)n]/mb. 

We see that in the FM phase the gauge fields acquire a 
"mass" term through the Anderson-Higgs mechanism. In 
other words, the excitations corresponding to d acquire 
a finite energy gap Eg ^ nia/kp, up to a multiplicative 
constant of order one, where kp is the Fermi momentum 
of holons. Using the mean-field parameters, m/ can be 



estimated as m^- 



(1 — x/4:)tl , and wc have 



En 



AV2Do/{tP) 



t ^3 + W/f{x) 



4 



1 - fi^) 



2Do, 

(45) 

where W = ep/t and f{x) = l — x/A. Equation (|45)) gives 
rise to the values Eg/t = 0.1273 — 0.2381 in the doping 
range 0.15 < a; < 0.5. For t ~ 0.3 eV, this corresponds 
to Eg ~ 0.038 — 0.071 eV, which is about the same order 

of Te. 

When the energy is much lower than Eg, one may fur- 
ther integrate out a, the resulting effective Lagrangian 
can be written a,s Cfm = Cq + Cmt where 

Co - h^idr-pf)h+^\Vh\^ + Pdrb + Do\Vb\'', (46) 



2m/ 



and 



|(j/.+j6)' + f (Vp, + Vp,)2 . 

In Eq. (j47p . we only keep the most relevant operators 
around the FM fixed point, described by Co, in the sense 
of RG, and J ft (jb) and Ph{pb) are the current and density 
operators for the h (6) fields, respectively. 

Since Cint is irrelevant around the fixed point described 
by Co, our low energy effective theory shows clearly 
that the FM phase of the QDE model exhibits the phe- 
nomenon of spin-charge separation. That is, the low en- 
ergy excitations are holons described by the h field, which 
carry charge e and spin zero, and the FM spin waves 
(magnons) described by the b field, which are charge neu- 
tral and carry spin-1. This is very different from the 
usual itinerant ferromagnets, where the elementary ex- 
citations are the dressed electrons, which carry charge 
— e and spin-1/2, and the magnons, which are collective 
excitations in the particle-hole channel. 



(47) 



SnibUd 



2mb 
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V. CONCLUSION AND DISCUSSIONS 

To summarize, we have studied the low-energy physics 
of the QDE model based on the slave fermion formula- 
tion. The most important feature of this approach is 
that the effects of large values of U and Jh are taken 
into account right from the beginning. (For the one- 
orbital model, the large U, Jh limit is implemented by 
the NDO condition.) This results in a non-Fermi-liquid 
ground state, in contrast to most of the previous stud- 
ies on the DE model. A direct consequence following 
it is that the electron spectral function exhibits a broad 
asymmetric peak away from the Fermi surface, in addi- 
tion to the quasiparticle peak. Both our results about 
the low energy magnetic properties and low-temperaturc 
thermodynamics are also consistent with experimental 
data in the FM metallic phase. On the other hand, our 
prediction for the optical conductivity at low energy is 
only qualitatively consistent with experiments. Quan- 
titatively, the Drude weight we obtained is larger than 
that actually observed. This indicates that to have a 
good quantitative description for the low-frequency opti- 
cal spectra of manganties, some extra ingredients beyond 
the one-orbital QDE model must be taken into account. 
However, we should emphasize that, for such a scenario 
to be valid, the inclusion of these new ingredients should 
not affect the low-temperature thermodynamics and low- 
energy magnetic properties predicted by the one-orbital 
QDE model in any drastic way because the one-orbital 
model already gives rise to reasonable results on them. 

We have also studied the role of gauge fluctuations 
and derived a low energy continuum effective theory. 
A previous study in terms of the slave-fermion gauge 
theory claimed that the gapless longitudinal gauge fluc- 
tuations play a dominant role on the electronic spec- 
tral propertiesj2£ However, the longitudinal component 
of gauge fields is gauge dependent and physics should 
be independent of the gauge choice. Moreover, both 
the time and longitudinal components of gauge fields are 
screened in the metallic phase, which cannot affect low 
energy physics dramatically. What we do in Sec. IIVI is 
to explicitly show that the gauge fluctuations in the FM 
phase, both the transverse and longitudinal components, 
are completely screened due to the Anderson-Higgs mech- 
anism. The resulting electron spectral function in the 



FM metallic phase consists of two parts: one the sharp 
quasiparticle peak at the Fermi energy and the other an 
asymmetric broad peak away from the Fermi energy, in- 
stead of a broad quasiparticle peak at the Fermi energy. 
These predictions can be observed from the ARPES or 
STS measurements. 

From the theoretical point of view, the slave fermion 
formulation also gives a clear picture of the apparently 
observed Fermi liquid behavior and its possible devia- 
tions. Within this framework, the Fermi-liquid- like be- 
haviors, such as a small value of the Drude weight in 
the optical spectra and a small coherent quasi-particle 
peak in the electron spectral function, are consequences 
of the Bose condensation of the spinon field. Therefore, 
the spectral weight is directly related to the magnitude 
of the condensate. Based on the similar reasoning, we ex- 
pect that in low dimensional or layered materials where 
the quantum fluctuations tend to kill the Bose conden- 
sate, the coherent quasi-particle peak will cither disap- 
pear completely or be greatly reduced in comparison with 
the three-dimensional results we obtained hcrc .^^'^^ This 
is actually observed in the ARPES measurements for lay- 
ered compoundsJ^ii^ Another minimal extension of the 
one-orbital QDE model, which has the potential to re- 
duce the quasi-particle peak or the Drude weight, is to 
include the orbital degrees of freedom. In this case, we 
have two slave boson fields: one describing the spin fluc- 
tuations and the other describing the orbital fluctuations. 
If only the spin slave field condenses at low temperature, 
the strong orbital fluctuations will also tend to reduce the 
coherent parts in the optical conductivity and the elec- 
tron spectral function. Researches along these directions 
are in progress. 
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